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The input impedance of a cylindrical dipole in a homogeneous anisotropic ionosphere is 
determined for arbitrary values of the medium parameters and arbitrary orientation of the 
dipole with respect to the earth's magnetic field. A sinusoidal current distribution is 
assumed, as well as a low value of dipole excitation, so that the field equations may be 
assumed to be linear. The Green's function is expressed as a Fourier integral in spherical 
propagation-constant coordinates. In performing a residues evaluation of the radial integral 
in this coordinate system, it is found necessary to deform the contour differently in different 
angular regions. By an appropriate rotation of the coordinates, the boundaries of these 
angular regions are made to depend on a single angular variable. An expansion of the result 
in power series in the small parameter e= (dipole radius/dipole half-length), in which only 
terms in log(l/e) and 6° are retained, leads to considerable simplification. It is found that 
the two dominant terms can be expressed as sine and cosine integrals, as in the free-space 
case, plus two finite single integrals which need numerical evaluation. The integrands of 
these latter integrals are expressed in terms of recursive routines. 

The numerical integrations are in terms of the polar angle of the wave propagation 
vector with respect to the earth's magnetic field. Singularities of the integrand are en- 
countered in the collisionless case when the applied frequency is below the plasma or gyro- 
frequencies. It is then necessary to carry out the numerical integration in the complex 
plane. Consequently the present numerical integration procedure is restricted to cases 
where this situation does not occur. Extensions of the treatment which remove these 
restrictions and methods of obtaining a more accurate current distribution arc; discussed 
briefly, the details being reserved for a later paper. 

Numerical calculations of impedance for selected values of the parameters are exhibited 
in a set of curves. 

1. Introduction 

For several years, ionosphere probing by rocket-borne probes has used the technique of 
the antenna impedance probe [Jackson and Kane, 1959]. The impedance probe has the special 
merit that it is affected by the ionospheric parameters in the immediate vicinity of the probe — 
roughly in the order of a wavelength or less of the probe excitation. Consequently, this type 
of probe gives local values of the parameters, rather than integrated values as in the Seddon 
[1953] type of experiment. It also can be operated at an excitation level low enough to avoid 
the creation of additional ionization in the medium by the rf field. 

The proper interpretation of the impedance measurement of an antenna probe requires a 
knowledge of the input impedance as a function of the medium properties. The type of antenna 
used in rocket probing can be rather closely approximated as a cylindrical dipole. Fortunately 
one is interested in the impedance change produced by the ionized medium, so that deviations 
from the ideal cylindrical form can be expected to play a secondary role. For this reason, 
the subject of this paper is the calculation of the input impedance of a cylindrical dipole in a 
magneto-ionic medium. 

Rocket experiments have shown that the antenna acquires a negative charge. This nega- 
tive charge repels the electrons immediately around the antenna, with the result that the antenna 
environment is rendered inhomogeneous. This effect is a very important one from the stand- 
point of the application of the impedance probe technique to the ionosphere. If this type of 
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probe is to be useful, the nature of the inhomogeneity produced, and its effect on the antenna 
impedance must be capable of determination. The former of these two problems appears to 
be the more formidable one at present. In any event, a solution for the homogeneous case is 
a necessary first step, both to develop an insight to the nature and magnitude of the effects 
produced by various values of the medium parameters as well as to serve as a basis for a possible 
perturbation technique for the inhomogeneous distribution. The case of a homogeneous 
medium will be assumed in this paper. 

The calculation of input impedance of a dipole is a straightforward problem when the 
distribution of current over the dipole is known. The current distribution is not arbitrary, 
however, since, in principle, it can be derived from Maxwell's equations and the known boundary 
conditions at the dipole surface. The key difficulty arises when one attempts to determine this 
current distribution from Maxwell's equations, since one is confronted with the problem of 
solving an integral equation [King, 1956]. Various iterative methods have been employed, but 
the accuracy of the result, as judged by comparison with experiment, is sensitive to the technique 
used. 

In first approximation the current distribution along the dipole is sinusoidal. This ap- 
proximation is quite good for very thin dipoles. The sinusoidal distribution may be considered 
to be the result of guided waves propagating in the medium along the outside of the conductor, 
at a velocity appropriate to the external medium, and perfectly reflected at the open ends, the 
interference between the two oppositely directed wave trains resulting in a standing wave with 
zero current at the outer ends. 

Fortunately, as already mentioned, in the case of a dipole used as an impedance probe we 
are not interested in the exact calculation of the dipole impedance, but rather in the impedance 
change upon entry into the ionosphere. This change is thus a difference quantity, so that small 
deviations in the nature of end effects which are occasioned by the use of only an approximate 
current distribution can be expected to largely cancel out when the difference, or change, from 
the free-space value is formed. Furthermore, the results obtained on the basis of a sinusoidal 
current distribution may be made the basis for a variational procedure whereby the current dis- 
tribution is optimized, so that the initial impedances values can be refined. In this paper we 
assume a sinusoidal current distribution, reserving the variational treatment for a subsequent 
communication . 

A further assumption is made that the amplitude of the motion of the free electrons in the 
medium in response to the electric field of the dipole is so small that the refractive index is given 
by the standard Appleton-Hartree formula. This assumption makes the field equations linear, 
so that Fourier resolutions are admissible. 

In the treatment given in this paper, the input impedance will be obtained by equating the 
complex power passing from the surface of the dipole into the medium to the complex power 
supplied to the dipole at its input terminals. In this formulation the Green's function is ex- 
pressed as a Fourier integral in spherical propagation-constant coordinates. In performing a 
residues evaluation of the radial integral in this coordinate system, it is found necessary to de- 
form the contour differently in different angular regions. By an appropriate rotation of the co- 
ordinates, the boundaries of these angular regions are made to depend on a single angular 
variable. An expansion of the result in power series in the small parameter e = (dipole radius/ 
dipole half-length), in which only terms in log (1/e) and e° are retained, leads to considerable 
simplification. It is found that the two dominant terms can be expressed as sine and cosine 
integrals, as in the free-space case, plus two finite single integrals which require numerical 
evaluation. The integrands of these latter integrals are expressed in terms of recursh e routines. 

The numerical integrations are in terms of the polar angle of the wave propagation vector 
with respect to the earth's magnetic field. Singular points of the integrand are encountered 
in the collisionless case when the applied frequency is below the plasma or gyrofrequencies. It 
is then necessary to carry out the numerical integration in the complex plane. Calculations 
given in this paper are restricted to cases where this situation does not occur. 2 The closed- 
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form terms, however, are valid in all parameter ranges. 

Numerical calculations of impedance for selected values of the parameters are exhibited in a 
set of curves. 

2. Formulation of the Problem 



In view of the preceding discussion, we now undertake the calculation of the input imped- 
ance of a cylindrical dipole in an infinite homogeneous ionosphere having a constant super- 
imposed magnetic field. This calculation is based on the following assumptions: 

(a) The current distribution along the (center-fed) dipole is a sinusoidal function of dis- 
tance along the dipole, being zero at the outer ends, and is uniform around the circumference 
of the dipole. 

(b) In the calculation of the field at a typical point in space, the current can be considered 
as concentrated along the axis of the dipole. 

(c) The ambient ionosphere is homogeneous, and its refractive indices are given by (he 
Appleton-Hartree relations. 

The dipole, of radius r and length 21, is considered to be fed at its center, with a current 
distribution as shown in figure 1. 

We shall find it convenient to employ two coordinate systems, 2 and 2'. Coordinate sys- 
tem 2 has its 2-axis along the earth's magnetic field, II . For 2', the s'-axis coincides with that 
of the dipole, and makes an angle 6< tt/2 with z. The relative orientation of 2 and 2' is chosen 
so that the s'-axis lies in the ?/z-plane. Thus the z-axis of 2 and the x'-axis of 2' coincide. The 
orientations of the two sets of axes are shown in figure 2. 

We shall assume a time dependence of e io>t for all field quantities and sources, national- 
ized mks units will be used throughout. 

The impedance will be determined by equating the complex power (i.e., volt-amperes) 
supplied to the antenna at its feed point to the complex power supplied by the dipole to the 
external medium, where 4 the latter is obtained by integrating the normal component of the 
complex Poynting vector over the surface of the dipole. The result is 



Z= 
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Figure 1. Illustrating dipole coordinates and current 
distribution. 



Figure 2. Coordinate systems. 

x,y,z is a rectangular system oriented with the earth's magnetic field, 

while x',y',z' is oriented with the dipole axis along z' '. 
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For the calculation of the integral in (1), we introduce cylindrical coordinates (r' , cp', z f ) 
in 2', where <p' is measured from the x'-axis. If we denote unit vectors by e, then 

ds=dAe r >. 

The actual current flows on the surface r'=r . However, for the calculation of the fields 
we make the usual assumption that the current is a line source on the axis of the cylinder, 
r' = 0, and directed along the axis. 3 Thus, if j(x') represents this line current, 

#tan K (r , ^ 2 0=^J(20X£ r ,=~rJ(z')V. (2) 

On inserting (2) into the integrand of (1), we obtain 

Zz= ~2^P\ fj j(z ' )(Ez C ° S 6+Ey Shl 6)[Wdz '' (3) 

The problem thus reduces to a determination of the electric field E. Once this has been 
found, the values of E y and E z can be inserted in (3), and this then evaluated to find the input 
impedance Z. 

3. Calculation of E 

We now turn to the problem of finding the electric field E. This obeys Maxwell's equa- 
tions for the medium. In view of the free electrons and the superimposed static magnetic 
field of the earth, the medium is characterized by a dielectric tensor, so that a given component 
of field is due to current components in all three coordinate directions. On the assumption 
that the field equations are linear (implying sufficiently small vibrations of the free charges), 
Maxwell's equations guarantee that a solution must exist for E(x) of the form 

E i (x)=jG ij (x\x,)J j (x 1 )(Px 1 (4) 

where cPx x is a compact notation for dx lx dx ly dxi z , and the convention of summing over repeated 
indices is used. 

In (4), GijixlXx) is the (i, j) component of the Green's tensor. Physically, it represents 
the electric field in the ith. direction at x due to the jth component of the source current density 
J at x lt Thus GijixlXi) propagates the effect of J at x 1 to the point x. 

A derivation of the Green's function GijQxXi) has been given by a number of authors (see, 
for example, Kogelnik [I960]). The result may be written in the form 

^■(x|x 1 ) = -^J^^ e - i - (5) 

From (4) and (5), 

E t (x) = -*j0 JVz, f<Pq ^ J,(xde-«-> (6) 

In (5) and (6), p is a numerical distance which is defined by 

p=k (x-x 1 ) (7) 

# being the free-space wave number. 

In (5) and (6), A Jt is the (j,i)ih cofactor and A is the determinant of the following matrix: 



M(q) = 



— q 2 (Li—K2i q 2 —q 2 2—K 22 —q 2 q 3 —K 23 



-qzqi—Kn ~ ^2— K32 q —qs—KsU 



(8) 



3 It can be shown that this leads to errors of order no higher than terms neglected in the present treatment. 
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in which 



(L 2 =g 2 i+q 2 2+g_l 



(9) 



The Kfj are components of a tensor, which, in the coordinate system 2, takes on the partic- 
ularly simple form 

(ot\ —ia 2 0\ 

ia 2 a, 1 (10) 

P aj 



/herel 



which are interrelated through 



1 xu 

Oii=l o 2 

a 2 = ^ 2 

u—y z 
1 x 



al=(ai—a 3 )(a L —l) 



(11) 



x, y, z arc the usual normalized (plasma frequency) 2 , gyrofrequency, and collision frequency, 
respectively: 



and 
Then 



5C=(w,v/o>) 2 ^ 

y=a) H /u: 
Z =v/a) 
u=l—iz j 

A(q) =det {M(q) } = -<*(ai-oi) (ql-al) 
where <x\ and a% are the roots of the biquadratic 

a 3 <r 4 -[2a 1 a 3 - (ai+« 3 )pV+ (p 2 -a 8 ) (a^-af +<£) =0 



with 



These roots thus are given by 

1 



01,2=2^- t2a ia3 — (ai + a 3 )p 2 ± [(a 3 — «i)y— 4a 3 a2p 2 +4aIai] 1/;; } 



(12) 

(13) 

(14) 
(15) 

(16) 



Since (3) does not involve E x , so that j x is not required in (4), we need only four of the nine 
matrix elements of Gtj(x\xi) for the present problem. These are G 22 , G 23 , G 32 , 6? 33 . Then from 
(4) and js=jV cos 0, jy—jz' sm #> (3) becomes 

+cos 2 ^33(x , |xn+sm^cos^c? 23 (x / |x; , ) + G r 32(x / [0]}. (17) 

4. Evaluation of Z 
To evaluate Z, we introduce dimensionless variables by putting 

(18) 
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where k a is the propagation constant associated with the current distribution. For generality, 
we assume that the propagation constant k a in the current distribution may be complex, so that 
A likewise may be complex. Then (17) becomes 

(19) 

Formally, the final step is to substitute the required values of 6? z; (x'|x") from (5) into (19) 
and perform the integrations to obtain the value for Z. The quantity in braces in (19) involves 
the ratio A r (q)/A(q) where A(q) is given by (13), and N(q) by 

2V(q)=sin 2 6[(q 2 — a,){(f 2 — a 3 ) + g?(«3-ai)] + cos 2 6>[( 5 2 -a 1 )(^-a 1 )-o:|]+2 sin cos Oq 2 q z {q 2 —a x ). 

(20) 
Consequently (19) becomes 

z -dyi!>wi>i>T* w >> m <-'''' <2,) 

in which f an d k , respectively, are the free-space impedance and propagation constant. 

We now are ready to introduce the assumption of a sinusoidal current distribution. We 
express this distribution, which is illustrated in figure 1, in the form 



j(f)=/„sin (A-f), 0<f<A^ 
=/ sin(A + f), 0>r>A J 



f- (22) 



and use this form for both j*({'i) and j(£' 2 ') in (17). As noted earlier, we assume that the 
propagation constant k a in the current distribution may be complex. 

With the assumed form of current distribution (22), it is possible to carry out all the 
inner integrations in (21). The quantity p in the exponent in (21) is the numerical distance 
between the current element on the axis at f" and the current element on the surface of the 
cylinder at f'. Hence 

q-p=q_AV — £")+q V 'R sin <p'+q x ,R cos /, 
where 

R=k r . (23) 

When the exponent in (21) becomes 

— i[<lz'({' — ti)+R(qx' cos y+qy sin if')]. 
By putting 



q x > cos v'+qy sin <p , = -y/q 2 x , + g 2 y , cos (<p'+y), 
y=t^n- 1 (q r /q x/ ), 
the innermost integral in (21) becomes, in view of the 27r-periodicity of the integrand 

„'+ v sm „o ^ =2 7rJ (i2V^+^)- (24) 

In the 2 system the Bessel function in (24) becomes 



f»27T 

Jo 



p — iR(q x , cos ^ 



Jo(R^ql + (q 2 cos 6— q 3 sin 0) 2 . 

Although (24) has been evaluated in closed form, it will be necessary later on to revert to the 
integral form in carrying out a residues integration. 
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The evaluation of the remaining integrals in (21) is straightforward, so that we obtain 

y '= r .(f.f^'£r *•>■*•>*"" (25) 



in which 

and 
In (27) 

in which 



C= 3 |' r ° 2A i (26) 
2ir I sin 2 A | 

^(q)=«/o(fiV2i + («2 cos 6-q z sin 0) 2 ) ■ GG*. (27) 

„_ A cos A (g 3 cos d J r q 2 sin 0)— cos A 

A (g 3 cos0+g 2 sin0) 2 -(A/A o ) 2 ( ' j 

A =W. (29) 



Study of (26) shows that a transformation from rectangular ^-coordinates to spherical 
coordinates is advantageous, since then the radicand in (18) does not pass through zero for 
real angles. The radial integration in this spherical coordinate system is then carried out by 
residues. One then gets analytically different results in different angular regions. In order to 
simplify the division of these regions, transformations (rotations) of the spherical coordinates 
are introduced in order that the regional boundaries depend on only one of the angular variables. 
The analytic forms of the integrand are then written out in full for each region, and further 
combinations of the component parts of the impedance expression can be made in order to 
simplify the subsequent manipulations. 

The integrands are then studied as power series in the parameter e=(dipole radius/dipole 
half-length), which is a small number for practical dipoles. The dominant terms of these 
series arc 4 examined as functions of the angular variables. It is found that the two dominant 
terms can be expressed as sine and cosine integrals, as in the free-space case, plus two single 
integrals which need numerical evaluation. The integrands of these latter integrals can be 
evaluated recursively. It is believed that integrals similar to these are inescapable because of 
the form of the propagation constants in the magneto-ionic medium, which contain the angular 
variable under a radical. 

In sections 6 and 7, the consequences of assumptions (a) and (b) above are examined 
briefly. With respect to (a), it is pointed out that variational procedures can be formulated 
which lead to a determination of a more precise current distribution. With respect to (b), it is 
reported that this assumption causes no error in the dominant terms in e, retained in the present 
calculation. 

Following the method of analysis outlined above, (25) is first transformed to spherical 
^-coordinates by introducing the following transformations: 



qi=q sin /z cos j8 
q 2 ~- q sin ii sin /3 
q 3 = q cos /z 

q = (ql + ql + ql) 1/2 



(30) 



(/<l\<lq><lq3= dp sin fidy q 2 dq. 



It is evident from (3), (4), (11) and (12) that the integrand is an even function of q. Hence 
we can write 



/» oo i r* oo 

Jo •••^<n>L---^- 



385 



This will facilitate an evaluation of the ^-integration by residues. 

In the subsequent residues evaluation of the ^-integration, it is necessary to expand the 
function F(q) in (25) into a form containing exponentials. The Bessel function in F(q) can be 
written as (see (24)) 



1 C 2t 



where 



<£— sin jjl cos (3 cos <p-(-(sin /* sin f$ cos 6— cos n sin 6) sin <p. 
Using the transformations in (30), (28) becomes 



G= fi 



cos (A rg) — cos (A Q) 



where 



tY-12 2 

r = cos jjl cos 0+sin /x sin sin /3, 
&=A/A . 
If the above changes are introduced into (25), it becomes 

Q (*2x f2* (*v r™ N 

Z =T~ d <P dfi sin ndfx q: 2 ~GG^e- iR ^dq. 

47rJo Jo Jo J -co A 



(31) 
(32) 

(33) 

(34) 
(35) 

(36) 



The functions <£ and 6767* control the convergence of the integrand at |#|->°°, and these 
depend on the angles <p, n, (3 through (28B and (30). In order to simplify this angular de- 
pendence, we now make rotations of these coordinates in order that the boundaries of the 
regions of convergence depend on only one angular variable. This is accomplished with the 
transformations 

cos ix cos 0+sin M sm p sin 0=cos a ^ 



sin n cos /3=sin a cos \p 

sin \x sin /3 cos 0— cos jz sin 0=sin a sin ^ ^ 

followed by the transformations: 

cos ju^sin a sin X 

sin yu' cos ft' =sin a cos X 

sin \x r sin /3'=cos a 



(37) 



(38) 



In this process 



dp sin ptrf/z changes to d(3' \ sin //rf/z'. 

JoJo Jo Jo 



Then (28) and (30) become, respectively, 

<£ = sin // cos /3' 
r = cos cr^sin /x' sin 0'. 



(39) 
(40) 



From (39) and (40), it is seen that $ and r have the same //-dependence, so that the convergence 
of the ^-integration will depend only on the angle jS'. 

The various coordinate transformations are depicted in figure 3. 

Upon introducing (30) into (13), A may be written as 



A=-[a l + (a z -a 1 )cos 2 f ,](q 2 -ri)(g 2 -7i) i 
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(41) 




X. LC 



Figure 3. Coordinate system showing the various angles used in transforming valuables. 

,y,z is a rectangular system oriented with the earth's magnetic field 1 1 ... and n,v,w is a rectangular system oriented with the dipole. Inset shows 

numerical distance p. 



where 



2= 2a 1 a 3 +(a l —a 3 ) sin 2 /xdb Jj 
2[o; 1 + (a 3 — ai)cOS 2 //] 



(42) 






xy 



u 2 —y 2 



r v 2 ~i l/2 

4a!cos 2 //+^2sinV 



(43) 



In (43), the x, y, u of (12) have been reintroduced for compression and in anticipation of 
numerical calculation. It should be noted that, for ^= 1(2=0), the radicand is positive and 
has no zeros for real /j,0</x<V/2. For future reference we write down 



^2 1 „2_ 2a 1 o£3+(tti— «3)sin 2 /i 

1 1\ ' 2 — 1/ \ ^^ 2 

a 3 (a 2 a 3 J r ol\—olz) 



rlrP 



ai + (a 3 — «i) COS 2 ju 
In a similar way, (20) for N becomes 



(44) 
(45) 

(46) 
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where 

Jf= cos 2 <t=t 2 (47a) 

{§ =«! cos 2 0-+ (ax— a 3 ) sin 2 sin 2 o- cos 2 \f/+ (ai cos 2 0+«3 sin 2 6) 

=a lT 2 + (ai-a 3 ) sin 2 cOS 2 ^(l - r 2 ) + ( ai cos 2 d+a s sin 2 0) (47b) 

C =«i«3+ («i— as) cos 2 (47c) 

with the above transformations, (36) takes the form 

Z= _^ r d + fV T **"'*»' r q2NG °* e-^da- (48) 

Z 4^J ^Jo ^Jo«i+(a3-«i)cosVJ-»(3 2 -rf)(2 2 -r 2 ) 2e " S ^ 

In the //, f$' variables, 

cos /*=cos a cos 0— sin <r sin sin ^ 



=sin// sin/3' cos#— Vl— sin 2 // sin 2 /3' sin 6 sin ^, (49) 

in which the radical is nonnegative, since 0<o-<7r. 

In order to perform a residues evaluation, it is necessary to split G and 6?* into exponentials. 
For compactness, we write 

%=R sin n' cos ft', 

rj=A sin // sin 0', 

rz=A/r), 

r 4 =-A*/v. (50) 

Then 

0G * e -iB*<Jl^ g(g) (51) 

r 4 ((T-r 2 3 )(q 2 — rl) 
where 

^(g) = J [ 6 -*fl^+^> 4. 6 -^<«-2,) _ KcosA+cosA^fe-^+^+e-^^-^l+Ci+cos Acos A*)e"^. 

(52) 
Then (48) becomes 

4ttJ Jo Jo ai+fe— «i) cos 2 m 

in which 

n Mf gW(g)ff(g)tf g 

V r 4 J_„ ( 3 2 -r?) 
The factor 



. (2 2 — ^)(g 2 — ri)(q 2 — ^)(2 2 — *1) 

ai+(«3— ai) cos 2 \l 



(54) 



in the denominator will be left in that form, since ultimately we shall return to the <p, n, /3 coordi- 
nates in making use of this factor. To facilitate the evaluation of the ^-integral of (54) by 
residues, we have written H(q) in (52) in exponential form. 

The poles ±r< of (54) occur in pairs; the contour passes below — r f and above +r f ; this 
choice makes e~ iuD/ an outgoing wave for distance Z?'>6, since the r t have negative imaginary 
parts for 2>0. Writing the integral so as to include all poles, the g-integration consequently 
yields 

where j, k, I with j^k^l, represent the three other subscripts possible after i is chosen. r t is 
chosen as that root of (42) which has negative imaginary part. H(r t ), in view of the discussion 
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above, is given by 

H(r t ) = l[e- ir itt +2l >\+e- ir ttt- 2l >l]-%(cos A+cos A*)[e-' r 'ls + *l+e-* r *'*-'i] 

+ (i+cos A cos A*)e~ ir iW. (56) 

The absolute value signs in the exponents require that // and fi' be real in the remaining 
integrations. 

The relative values of £ and rj may be divided into three regions, in which the integration 
over (3' must be treated differently. We define these by 

Region I: |f|>|2i,|, (57a) 

Region II: \2rj\>\^\>\rj\, (57b) 

Region III: |r?|>|£|. (57c) 

From (56), we may write after some rearrangement: 

in region I, Il(r i )=A(r i ) = e~ iT i l ^(cosr]r i —cos A)(cosr]r i —cos A*) (58a) 

in region II, H{r l )^A{r l ) + \[e- ir i\ 2 ^ cos (/%•£) — e~ ir ^\ cos (2r^)] 

=A(r t )~ sin [r^rjl-m 

^Aird+Biu) (58b) 

in region III, H(r i )=-4(r / )+S(r,)+i(cos A+cos A*) sin[r,(|i|| — 1$|)] 

^(r0+5(r<)+(7(r,). (58c) 

Symbolizing the angular integration over region I by I , etc., and the integrand by S(r t )H(r t ), 
we have 

Z= f A(r i )S(r i )d(l'+ f B(r i )S(r i )W+ f f7(r<)S(r,)<i/5' 

Jl + H + III .7 II + IH Jill 

Z, + Z 2 + Z 3 (59) 

The boundaries of the regions I, II, III, which depend on fi' , will be discussed in section 4.1. 

4.1. Evaluation of Z\ 

Since the regional sum I + 11+ III in Z x comprises the full range of angular variables of 
integration, we may conveniently rotate back to the original <p, /*, (3 variables in (36) to take 
advantage of the fact that the complication of the radical in (43) involves only ^. Also we 

4 

note here that 4(r 7 )S(^) has no singularities for ^=3, 4, so that 2 m (55) may be replaced 

2 

by T] m evaluating Z Y . Hence we obtain 

'\=\ 



J» 2ir (* 2ir (* T 

<h dp 
o Jo Jo 



sin ijidfji 



a! + (a 3 — <xi) COS 2 \x 



1 f at-/ \ /cos '//'i— cos A \ /cos nrj— cos A*\ _, r . 



■r,N(r a ) 
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/cos rjr 2 —c 



ryr 2 - cos A \ / cos ryr 2 -cos A* \ ir M \ , . 
7? 2 -A 2 A ^ 2 -A* 2 / / K } 



in which 

7T^|A 2 |C 



Ci=- 



47T 

A 



= -60(//X ) 



sin A 



(61) 



Naturally, this form for Z x is obtainable from (36) by treating the two factors G, G* as entire 
functions whose growth far from the origin and from the poles is dominated by e~ ig * for |£|>0. 

Now the roots r h r 2 of A are numbers of order unity, and are exactly unity in the free- 
space case. The maximum value of |£| is R, the dipole radius measured in free-space wave- 
lengths, which is a small number for thin dipoles. In that case, then, we can expand in powers 
of R and expect rapid convergence. The leading term in this expansion is obtained by setting 
e~ ir i.2W = l. For this leading term, then, <p vanishes from the integrand, so that the ^-integra- 
tion yields a factor 2t. 

Since there are functions of /x buried in the forms for r 1} r 2 , no analytic /x-integration is 
in sight. Hence we resort to numerical integration for this, restricting /x to real values. 
Since the r t do not involve /?, only the /^-integration needs to be done numerically. 

In order to prevent the ( ) factors in (60) from assuming the form - at some point in 

the p, jit-range, we expand these factors in series. A typical factor may be written as 

cos a— cos A _ ~ (— l) n a 2n — A 2n 
a 2 - A 2 ~hi (2n)\ a 2 - A 2 ' 

Then, after some manipulation, we obtain 

/ cos q-cos A \ / cos q-cos A* \ » 2 

V a 2 - A 2 A a 2 - A* 2 )~hi> n) ( } 

D "=(~V n g § § [2 ( m +j+l)]\[2(n-m+k+l)}\ (63) 



where 



D n , which involves only the constants A, A*, can be calculated and tabulated against n in 
advance of the main calculation. Since 

a = r i r} = A r i cos (r=A 7 , i (cos jjl cos 0+sin n sin /? sin 6), 

the a 2n lead to terms containing integral powers of cos 2 a. Furthermore, as can be seen from 
(40) and (46), N{r t ) involves a and (3 in the forms cos 2 a and cos 2 f3. Hence in the /^-integration 
we need to evaluate 

»2ir 

cos 2n adfi, (64) 



and 



2tt Jo 

i r 2ir 

Qn=2 T cos2 £ cos2n «#■ (65) 



These may be evaluated to obtain 



P »=S (2m)![ ( 2?r!-m)] ! (c0S " C0S ^ ( "" (sin * sin ^ < 66 > 

«•=& (2m)![2(,-m)] ! (c0S " C0S ^"^ (sin " sin ^ < 67 > 
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where 

(2m)! 



Xf n o2«*y,l«y. ?' y* 



m 2 2m m!m! *"" 2(m + l) 

2 n 1 

(x n may be calculated recursively from £o=l> ^« — ~~^ — #»-i-) 

With these results, we get from (60) on putting 

e - ir iW=l+0(R) 
and using (46) for N(r t ), 

Jo Jo ai+(a 3 — ai)cos> fifTi r^—rj 7 fr 

Jo Jo «i+(tt 3 — ai) cos M 

- jy (r? » + 5_ rr? + 5) _ @ (rf + 3_ r 2n + 3 )+ g(r 2, + l_ r 2, + l ) 

• 2j — -2 — ^ — T *v Ov 

71 = /l I 2 



The sum mat ion contains terms of the form 

^J '2 ^n 

ri—rl ~r l -\-r 2 
where 



(69) 



c„=(r?+r|)c w _ 1 — r?r|c re _ 2 . (69a) 

Hence, putting cos /*=£, we obtain 
eft 1 



Jo 



oLi + (oL Z —ai)t 2 r\+r 2 

00 

• S ^AS tt {fc + 2-«A+i)P w+ i + [(«3-«i) cos 2 0— a 3 ]c»xiAi 

7J = 

+[a 1 a 3 - (az-ca) cos 2 W^ + fe-^K+i sin 2 0(l-* 2 )£ n }. (70) 

It should be noted that the appearance of 6, both explicitly and implicitly in P n , Q tl , is always 
as an even power of cos 6 or sin 0, so that 

Z 1 (6) = Z l ( 7 r-d), 

as is required physically. Equations (42), (63), (66), (67) and (69) permit the recursive evalua- 
tion of the integrand of (70) for each /*. Thus what remains is a numerical integration over /x. 

4.2. Evaluation of Z 2 

In Z 2 (and later in Z 3 ), we leave \f/, y! ', /?' as the angular integration variables in order 
to have the boundaries of regions II and III determined by the single azimuthal variable /?'. 
Later, after introduction of suitable approximations, it will be possible to return to the original 
variables <p, fx, /?. 

The integral Z 2 is taken over the region 11 + III, which, from (57b), is given by 

m<i2ni. 

By defining 

_ 11 _r _l , 7t v 

e, -2A~2~r2 e ' (71j 
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this can be written as 

|tan/3'|>€. (72) 

The result of the ^-integration in (54) is given by (55). For Z 3 , we have to insert for 
H(r t ), in virtue of (58b, c). 

ff(r < )=J5(r,) = — ^ sin [r,(l2i, I — |f|)]: ^ sin [2r,|i|| (1— dl cot iS^I)]". (73) 

It contributes to simplicity to write Z 2 in terms of the quantities 

Si=ri\cos <r\=ri\sm p! sin $'\. (74) 



Then we have from (50) 

S 3 =A/A =ft 



S 4 =-A*/Ao=-Q*J 



(75) 



Si and S 2 , in view of the dependence of r x and r 2 on *z (see (42)), and the explicit factor [sin fi'\ 
in (74), are functions of 0', while S 3 and £ 4 are not. Then from (53), (55), and (58b, c), we 
have 



where 



1 4 9f S A (fi S' 2 (P 

G\= J_7 7~7~'Ij7o2 2\/o2 02W Q2 02\ ^ Sm [^ A Sj ( 1 ~ «1 COt /3') ] , (77) 



in which, with r=cos c again, 



® = (B =a l T*+(a-a i )sm 2 ecos 2 4 / (l-T i ) + (a l cos 2 8+a 3 sm 2 6) f" (78) 

C =T 2 <? =r 2 [a l a 3 +(a l -a 3 ) COS 2 0] 



r 



The term |cot fi'\ in the argument of sin [ ] in (77) and the |sin $'\ in the denominator of 
(76) are both symmetric about each of the principal axes in the /3'-plane. As function of S iy 
the summand in (77) expands in even powers of St, since S t sin [2A S f (l— €i|cot /3 f \)] is an 
even function of S t and the remaining terms are explicit in S 2 . Thus, in the summand, the 
cos a of (74) appears squared, and for this reason is symmetric about the principal axes in 

the /3-plane. The same cos 2 a symmetry appears in the factors J?] (B, G as given in (78). There 
remains the dependence of cos /x on £'; this appears both in the explicit factor \ I[ol\-\- {ol$~ a x ) 
cos 2 /x] and implicitly in Si, S 2 through their dependence on r h r 2 , which in turn depend on ji 
through cos 2 /jl, as is seen in (42). Now we may regard the integrand of (76) as expanded in 
powers of cos \i or cos 2 /x; one sees from (49) that such an expansion will involve odd powers 
of cos (j= sin \x' sin (3' , but that such odd powers of cos a will be multiplied in turn by factors 
odd in sin \p. Now \p appears otherwise only in the cos 2 yj/ term in (B ; thus we may visualize 
doing the ^-integral first, and we see that all terms in the expansion of the integrand having 
odd powers of sin \p give zero contribution to Z 2 . But these terms are exactly those which are 
odd in cos o-=sin // sin (3\ Thus, in view of the ^-integration, the terms in sin /?' due to im- 
plicit and explicit dependence on cos a are even, and therefore symmetric about each principal 
axis in the /3'-plane. Thus the entire integrand has this symmetry and we may henceforth 
confine /3' to the first quadrant and drop the absolute value signs. Then in (76), 

f . . .dp' =4 f W2 ...dp'. (79) 

Jii+iii Jtan-lei 
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The foregoing argument has also shown that part of the integrand of (76) which depends 
on the Si and on cos ^, either explicitly or via r u r 2 , can be expanded in power series in cos a= 
sin \i f sin f$' and that this is a series in even powers of cos a, in view of the ^-integration. Let 
us now regard this expansion as having been made. The ^'-dependence is then expressed in 
a series giving the ^'-integrals in the typical form 

fW2 dB' 

cos'mi-*! cot n m -£rR>' (80) 

J tan l *i Sill p 

Expanding (1 — e t cot (3 f ) m to get a double series, we have the typical integral 

h,r= f W2 (sin ^) 2n (e l cot pydfi'/mn 0'. (81) 

Jtan-'ei 

Now we take advantage of the fact that we are dealing with a "thin" dipole, one in which 
e = R/A is a small quantity. We wish to calculate terms of order log (1 /ei) and e?, and to neglect 
terms of order e? and e? log (l/e^ for n>l. By evaluating integrals of the type (81) with 
these considerations in mind, it may be shown that terms of the retained order occur only 
when ?i = 0, or r=(), or both, the last alternative giving rise to a log (l/ej term. We there- 
fore write 

Z 2 =Z' 2 +Z' 2 ' +0[d log (1/d)], 

where Z' 2 corresponds to n = 0, and /J/ to r=0, 717*0. Thus Z2 lias the integrand obtained 
from (76) by putting cos <r=0 in S t , in cos m, and in SC, (B,G , but retaining e, cot #' finite. 

From (78), we see that G is of order cos 2 a, and from (74) that SI and S'i are of this order. 
From (76), however, S 3 and S 4 are independent of cos a. Thus in the 2 of (77), for i=l, 2 the 
numerator is of fourth order in cos a and the denominator is of second order, so that the sum- 
mands vanish for i=l, 2 on setting cos o- = (). Thus we obtain 



2 'Jo Jo ai + (a 3 — «i) cos J M(. J».mii-'«i sin/3' 



^ [(X51- ®rS?)] sin [2A f S f (l-«, cot 0') ] ,o 9 > 

where, from (50), 

cos 2 /u =sin 2 sin 2 i/\ (83) 

and, from (78), 



Jxo— 1 

(B =(az— ai) sin 2 sin 2 ^+ai = a 1 + (a 3 — en) cos 2 ^ J 



i 



(84) 



We note that the integrand in (82) is free of //, and that the ^-dependency, in virtue of 
(75), (83), and (84), is of the form 



r* # = r 

Jo Q:i + (a 3 — on) sin 2 6 sin 2 iA J 



!7r d^ 

«i cos 2 \//+ (ai cos 2 ^+0:3 sin 2 6) sin 2 ^ 




=&r[a l (a 1 cos 2 0+a 3 sin 2 fl)]-**, 

from Bierens de Haan [1957] table 47, No. 10. After integration over \p and /*', (82) therefore 
becomes 

2-*Wir, A- £ '-<- c ° 8 ' '+5 y '-^-' ,in Ma-, eo, n. (85) 

Jtan -1 *! Sill p 1 = 3 AJiV^i — iOj) 

Writing 

0=€iCOtjS', 
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we now expand 

sin [2A #*(1 — 0)]=sin (2A Q S t ) cos (2A S^)-cos (2A S t ) sin (2A S t g) 

and further expand cos (2A Sig) and sin (2A Sig) in series. If we retain only terms of order 
log (1/ei) and €?, as before, we may sum the resulting series to give the result 

7>_ 2 pin2i ^ [«i(«i cos 2 ^+a 3 sin 2 ^)]-^?-l 

Zj 2 — TT 0|M | 2-J O / ■ Q2 02\ 

•/sin (2A S,) [tt(2A ^)+log ^^.]- cos (2A S«)Si(2A S r i ) \, (86) 

where Si (a;) and 6 y i(£) are the sine and cosine integrals [Jahnke and Emde, 1945], respectively, 
and in which use has been made of (71). This gives a closed form expression for Z' 2 . 

On the other hand, 2%. has the integrand obtained from (76) by setting eico\>fi' = in the 
final sin [ ] factor, so that 



where 



zj'=4c|a*| P# [* fa' r /2 -- ( — ^4 — — -^7> (87) 

2 ' ' Jo Jo Jtan-i«! ai+(<*3— «l) COS 2 fX sin /3' 

r _^ (^St- @ S 2 + (?)£, sin (2A S Z -) 

""0 >C_J /Qf2 Q2WQ2 Q2WQ2 02\ \P&) 



00 now is to be expanded in a series in powers of cos 2 <r, starting with the first power (since the 
term in (cos 2 a) already has been included in Z' 2 ). Thus in Z 2 the integrand has the leading 
term cos 2 <7=sin 2 // sin 2 /?'. This and all higher terms vanish like sin 2 ft' or faster near ft' = so 
that, even with the 1/sin ft' factor, the Z 2 integrand vanishes at ft' =0. This in turn permits 
the replacement of tan -1 e x by as the lower limit of the /^'-integral, since this introduces a 
negligible error of order e?. Hence (87) becomes 

Z' 2 '=\c\tf\ f 2T # f r <&i'f* , , G ° v — -^S>- (89) 

2 ' ' Jo Jo Jo ai+(a 3 — «i) cos 2 /x sin jS' 

After expansion of G in powers of cos 2 <7=sin 2 // sin 2 ft', we can make use of (38) and (37) to 
rotate back to the original coordinates <p, \x ) ft. This is desirable because then r l9 r 2 depend 
only on /x. The result is 

Z ' 2 '=i cm r d? f ' * ™*fr 2 . oo) 

4 Jo Jo |COSo-|a 1 + (a:3— a x ) COS 2 ^ 

^o in (89) is a symmetric function of the four Si. By expanding sin (2A #i) in a series, 
#0 becomes 

r A (^Sj-®ff+C)fr sin (2A fr) ^ ^ (-l^Ap) 2 ^ 1 ^Sj-gfct+Cff . B 

^-^ (^-^(^-^(SJ-flJ) ^ 2^ (2n+ l)! ( /S 2_ S 2 )(S 2__ S 2 )(S 2_ S 2 ) &< • 

(91) 
We can now make use of the results of symmetric function theory to obtain the desired expan- 

/»2ir 

sion of G in terms of cos 2 a. When this is done, the /3-integration |cos a \ 2n+1 dft can be devel- 
oped recursively, leaving only a /x-integration to be performed numerically. 
The expansion of a typical term of (91) 

§ (S?-SJ) (S]- SI) (S?-S?) = ^ 2 »+^ 
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is expressed recursively in terms of the coefficients A 2n of appendix A as follows: 

A =A 2 =A,=0, A Q =1, 1 

A2 n = (l\A2n-2 #2^2;i--4l*4-"-2n-6 ^4^2«-8 J 

di=^2 S'i=p +Pi cos 2 or 

d 2 =S 2' S r fS?=p 2 +PoPi cos 2 0- + P3 cos 4 a- 

d*=i: 2' S" S5SJS|=p 1 p a cos 2 (r+p P3 cos 4 cr 



rf 4 = II S 2 i = p 2 P3 COS 4 <r 



where 



(92) 



(93) 



(94) 



p = ft 2 +12* 2 , p 2 -12 2 l]* 2 =|l2 4 | ? 

P!=r?+r|, Pz=rtii. 

It is evident that the ^4 2 « will involve cos 2 <r, in general. We therefore write 

^2n=2 tt w ,„ COS 2 ™ a. 

71 = 

Hence we have from (92) and (93) 

a m ,o=a m , ,=a mt2 =0, for all m 
a 0t 3 = 1 , a m , 3 = for m 7* 

ttm,n = Po a 'm,n + \ P2flm,n-2~r PlQm-l,n-l PoPl^m-1, »-2 

H~PjP2^/«-l. w-3 P3 a m-2, n-2 T PoPz a m-2, »-3 P2Pz^m-2, n-i 

Then (91), together with (78), becomes 

£o=S S (-1)" £ A °! 2 il! COS^cr- {^.n+s-fl^n+gfa COS 2 g+a 3 sill 2 g) 
TO =i n=o (zn-f-i)i 

— a m>n+2 (ai— a 3 ) sin 2 sin 2 0- cos 2 ^— a w _ 1>w+2 «i+^-i, w +ikia 3 +(a:i— a 3 ) cos 2 0]} . (96) 

In (96), the factor sin 2 a cos 2 ^ has been retained in order to facilitate the rotation back to the 
original (<p, n, (3) coordinates. By the second equation of (37), this factor then becomes 
sin 2 n cos 2 13. Hence, upon rotation of coordinates, (89) becomes, in virtue of (90), 



(95) 



*'-Ti°ir*r 



sin \xd\i 



00 00 



«i+ (0:3—011) COS 2 \i 

(2A Q ) 2 ^ 



•S S i-^ T^hV} \™s°\ 2m ~ K {^+3-^+2(01 cos 2 0+a 3 sin 2 0) 

— dm, ^+2(0] — a 3 ) sin 2 6 sin 2 ju cos 2 $—a m -\, *+aai+a»-i. n+i[ai«3+ («i — o 3 ) cos 2 0] ] 

Since from (34) 

cos o-= r =cos /x cos 0+sin /x sin sin jS, 

the /3-integration involves integrals of the form 

Jn=~ f a '|a+6sini8|»rfi8, 

^7T Jo 

9n=<r f 27r cos 2 i8|a+6sin/3|V/3. 



(97) 
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These are evaluated in appendix B, again in terms of a recursive routine. Making use of the 
results obtained there, and making the substitution cos ju— t, (97) becomes 

Z' 2 '=t 2 C\Q 2 \ ,, ,. a >s S (-i) B 9 °L, 

'Jo «i + (a 3 — ajr ^o^=i (2^+1) I 

' { [« Wf n + 3 — «m, n+2(«l COS 2 + a 3 SHI 2 0) — a m _i, n+2 «l + «m-l, n+l(<*ia 3 

+ (ai— a 3 ) cos 2 <9)]/ 2 m-i+«m,n+2(« 3 — «i) sin 2 0(1 — t 2 )g 2m -i } . (98) 

The final step in the evaluation of Z' 2 ' thus is a numerical integration over t. 

4.3. Evaluation of Z 3 

The evaluation of Z 3 follows a procedure entirely parallel to that for Z 2 , it being only 
necessary to use 

e = R /A 

instead of e u and to use the appropriate value of H(r t ) from (59). It is easily seen that the 
only changes required from the procedure of section 4.2 is to replace 2A by A everywhere, 
and to multiply the result by the coefficient 

— 2(cos A+cos A*). 
Hence, writing 

Z 3 =Z 3 +Z 3 , 

we obtain from (86) and (98), respectively, 

7' o m\ r»2i / a i a *\ ^ [ a i(«i cos2 0+<*3 sin 2 6)]- 1/2 S 2 i — 1 
Z 3 =-27r 2 6 y | ti 2 \ (cos A + cos A*) J2 — ~ 07™ — ™^ 

•/sin (A SO [^(Ao^+log^-^-J-cos (AoS,)ff f (AoS«) V (99) 



Zr = -2T 2 C|fi 2 |(cosA + cosA*) f— 77-^ ^•SSH)" 

Jo «i + («3 — aijt «=om=i 



AJ 



2n + l 



(271+1)1 

• {[a m , n+z — a m , n+2 (ai cos 2 0+a 3 sin 2 0)— a m _ 1>n+2 c*i 
+a m - ltn+l (a l a s —(a z —a 1 ) cos 2 0)]/2m-i+am. w .+2(<* 3 — «i) sin 2 0(1 — £ 2 )# 2 m-i} . (100) 

5. Limitations of the Numerical Integration 

The numerical integrations to be performed in (86), (98), (99), and (100) in the variable 
£ = cos ju are carried along the real £-axis over the range <t < 1. From figure 3, it is seen that /x 
is the angle of propagation of the plane wave with respect to the earth's magnetic field. In the 
integrand, the term 

D= ai + (a 3 — a,)t 2 

in the denominator vanishes, when u=l, at some value in the range 0<£<1 when 

<2l#3^0. 

This is equivalent to 

P=(l-;r)(l-^ 2 )(l-z-2/ 2 )<0. (101) 

When D=0, the present numerical procedure does not give meaningful results. This is 
due to the fact that when the quantity A of (41), which is the denominator of the Green's 
function in (36), vanishes, so that the Green's function becomes singular at some point of the 
integration path. This difficulty does not occur when the collision parameter 2 is substantial, 
since then D vanishes well off the real £-axis. For small 2, or for P positive and small, a finer 
integration grid will also give valid results. 

396 



In order to avoid the singularity and obtain a single numerical process valid in all parameter 
ranges, it is necessary to carry out the numerical integration essentially in the complex £-plane.* 
This will be postponed bo a later treatment. 

The relations for the impedance components developed in this report, namely (60), (86), 
(98), (99), and (100), have been programmed for the IBM 7090 computer. The program has 
been satisfactorily checked against the free-space impedance in the high-frequency limit. 
Numerical calculations have been made for certain parameter values. These will be discussed 
in section 8. 

6. Propagation Constant k a and the Current Distribution 

As mentioned in section 4, the propagation constant of the assumed sinusoidal distribution 
of current along the dipole has been taken to be complex. This is reflected in the parameters 

A=kJ 
and 

a=A/A , 

which, consequently, are complex, in general. 

Our analysis tells us nothing about k a , since we have assumed the form of the current distri- 
bution, which amounts to assuming that we know k a . For trial numerical calculations, we 
have chosen values for k a near the average (over propagation direction) of the ordinary and 
extraordinary wave propagation constants. The value actually used was k a =^l—xk . 

In examining procedures for an optimum choice of k a , we have been led to a variational 
formulation of the thin dipole impedance problem. This type of procedure, introduced by 
Storer [1950] in the free-space case, leads to a method for finding optimum values for the current 
distribution, which is then used to obtain an even more accurate impedance estimate. Further- 
more, no major changes in the present analytical or numerical methods are required. These 
developments will be given in a later paper. 

7. Further Developments 

Our formulation is based on a Fourier integral representation of the anisotropic Green's 
function. Considerable effort has been expended, both by us and others, toward obtainingja 
closed-form representation of this Green's function. We have been able to show, however, 
that such attempts are in vain, and that a representation, such as the Fourier integral used here, 
is a necessity. 

Because we use a Fourier representation for the Green's function rather than a closed form 
as in the free-space impedance calculation, we have been able to evaluate precisely the error 
introduced by assumption (b) of section 4, namely, that the current can be considered as con- 
centrated along the axis of the dipole in the calculation of the field at a typical point in space. 
It turns out that the corrections for this assumption occur only in the neglected terms of order e 
or e log (1/e) or higher, so that this assumption causes no error in the dominant terms retained 
in this paper. 

Details of these developments will be given in a subsequent paper. 

8. Numerical Calculations 

In order to display the magnitude of the impedance as a function of the parameters, com- 
putations of Z were made for several frequencies. These were calculated for frequencies of 10, 
5, and 2 Mc/s and a set of ionospheric and antenna parameters supplied by Dr. R. G. Stone of 
the National Aeronautics and Space Administration. Table 1 lists the ionospheric parameters 
at the altitudes of interest. The corresponding values of x, y, z are plotted in figures 4 and 5 for 
10 and 5 Mc/s, respectively. From these it is seen that in an altitude region around 300 km, 
and of width increasing with decreasing frequency the parameter x becomes comparable with 

*Note added in proof: 
Suitable paths in the complex /-plane have been found for the numerical integrations Z x , z'/ , and 7.'/ . With a lG-point Gaussian numerical inte- 
gration procedure, calculations have been carried out at frequencies as low as about 2000 c/s without difficulty. 
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Table 1 



Altitude 


Electron 


Mag. field 


Collision 




density 




frequency 


km 


per cm* 


Gauss 


sec- 1 


100 


1.2X105 


0.52 


1000 


150 


1.7 


.51 


100 


200 


2.5 


.49 


500 


250 


3.5 


.48 


1000 


300 


5.2 


.47 


920 


350 


4.3 


.46 


840 


400 


2.7 


.45 


750 


450 


1.8 


.44 


670 


500 


1.3 


.43 


590 


550 


8.3X10 4 


.42 


510 


600 


5.7 


.41 


430 


650 


4.0 


.40 


350 


700 


2.9 


.39 


260 


750 


2.2 


.39 


180 


800 


1.7 


.38 


100 


850 


1.3 


.37 


95 


900 


1.1 


.36 


80 


950 


8.8X103 


.36 


75 


1000 


7.4 


.35 


70 
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HEIGHT, KILOMETERS 

Figure 4. Normalized ionosphere parameter versus 
altitude (based on data in table 1) ; frequency = 10 

Mc/s. 



/ \ x 

\ ~7^ — \__ Y 

1 / *\ \ " — " 
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HEIGHT, KILOMETERS 



Figure 5. Normalized ionosphere parameter versus 
altitude (based on data in table 1); frequency = 5 
Mc/s. 



unity or larger, so that the product P in (101) is small or negative. For smaller and greater 
altitudes than this region, computations of Z were made on the NASA 7090 computer. The 
resulting impedance, plotted against 0, the angle of the dipole's axis with respect to the earth's 
magnetic field, is plotted for a number of altitudes in figures 6 and 7 for 10 Mc/s, in figures 8, 9, 
and 10 for 5 Mc/s, and for an altitude of 1000 km in figure 11 for 2 Mc/s. 

The antenna in question is a half -wave dipole of radius-to-half -length ratio e = 1/300. Thus 
A=k a l has the value 1.5708 (i.e., tt/2). This corresponds to a situation where the length of the 
dipole is varied with altitude so that it would always have near-zero reactance. That is, the 
impedances calculated in figures 6 to 11 are based on the assumption of a sinusoidal current 
distribution of the form sin k a (l— \z\)=cos (k a z), k a being a real constant having the value 
-yjl—x k , independent of 0, where the ionospheric parameter x is that corresponding to the partic- 
ular altitude and frequency. 



398 



100 



90 



5 
X 
O 



300 km 



^350 km 



<-250km 
^400 km 



v 200km 
^450 km 



H 50, 100 krr 



,550km 600km 



r 



^500 km 



^-7nn ^ 



700 v 800 ^900 MOOOkm 



0° 15° 30° 45° 60° 75° 90° 



52 



48 



40 









\ 




300 km 


-~~-~~~ 


















— 
























































350 
250 




T__ 








200 
450 
500 


— ; — " 




___]___ 




1 1 




1 


■=■ 


550 
700 
1000 km 


t 












15° 










6 








)• 


30° 


45° 
B 


0" 


75° 90° 



Figure 6. Resistance of half -wave dipole versus 

angle to magnetic field; frequency = 10 Mc/s. 

Ionosphere parameters as in table 1. 



Figure 7. Reactance of half-wave dipole versus angle 

to magnetic field; frequency = 10 Mc/s. 

Ionosphere parameters as in table 1. 
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Figure 8. Resistance of half-wave dipole versus 

angle to magnetic field; frequency = 5 Mc/s. 

Ionosphere parameters as in table 1. 



Figure 9. Reactance of half-wave dipole versus angle 
to magnetic field; frequency = 5 Mc/s. 

Ionosphere parameters as in table 1. 
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Figure 10. Reactance of half -wave dipole versus 

angle to magnetic field; frequency = 5 Mc/s. 

Ionosphere parameters as in table 1. 
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Figure 11. Impedance of half-wave dipole versus 
angle to magnetic field; frequency = 2 Mc/s. 

Ionosphere parameters as in table 1. 
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Figure 12. Dipole resistance versus A = k a l; fre- 
quency =10 Mc/s, = 45°. 

Ionosphere parameters as in table 1. 
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Figure 13. Dipole reactance versus A = k a l; fre- 
quency =10 Mc/s, 6 = 45°. 

Ionosphere parameters as in table 1. 
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Figure 14. Dipole resistance versus A = k a l; fre- 
quency =5 Mc/s, altitude = 200 km. 
Ionosphere parameters as in table 1. 
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Figure 15. Dipole reactance versus A = k a l; fre- 
quency =5 Mc/s, altitude = 200 km. 
Ionosphere parameters as in table 1. 



In order to explore the magnitude of the change of impedance with physical length of the 
dipole, curves were also computed for the reactance and resistance as function of dipole length, 
for values of A = kJ near tt/2. The results are plotted here in figures 12 to 15, for the special 
conditions listed in the legends. 

In summarizing these plots, we first note that for 10 Mc/s and 1000 km altitude, where the 
ionospheric conditions have small effect and free-space conditions are being approached, the 
calculated impedance shown in figures 6 and 7 shows a negligible ^-dependence about a value 
identical with that of a half-wave dipole in free space [Carter, 1932]. 

Next, it is obvious from the graphs that the impedance is markedly affected in the altitude 
range where the ionospheric parameter x approaches unity. The impedance then becomes 
markedly dependent on angle, more especially the reactive component. This is especially 
noticeable in figure 10 for an altitude of 400 km. When this curve was first obtained, it was 
felt that these large values may have been the result of the approach of the parameter P of 
(101) to a small value (about 0.01) so that the numerical integration may have become unreli- 
able. To check this, the 16-point numerical integration was changed to 32-point, but the com- 
puted results changed negligibly. This is highly encouraging as to the accuracy of the numerical 
integration process . 

It should be pointed out that these curves are based on the assumptions stated previously, 
and therefore should not be taken as the best obtainable under the present theory. In par- 
ticular, the & a -value used here is real and independent of orientation angle, whereas the varia- 
tional estimates of k a referred to at the end of section 6 would lead to a complex ^-dependent 
value. But this refinement will be handled in a later treatment. 

It should also be pointed out that the present theory assumes a homogeneous medium 
around the dipole. It therefore does not incorporate the effect of the plasma sheath which 
may be formed around the dipole in an actual case. 

9. Summary 

In this paper, we have formulated the problem of the impedance of a thin center-driven 
dipole in a magneto-ionic medium by expressing the Green's function as a Fourier integral in 
spherical propagation-constant coordinates. In performing a residues evaluation of the 
radial integral in this coordinate system, it is found necessary to deform the contour differently 
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in different angular regions. By an appropriate rotation of the coordinates, the boundaries 
of these angular regions are made to depend on a single angular variable. 

An expansion of the result in power series in the small parameter e=(dipole radius/dipole 
half-length), in which only terms in log (1/e) and e° are retained, leads to considerable simplifi- 
cation. It is found that the two dominant terms can be expressed as sine and cosine integrals, 
as in the free-space case, plus two finite single integrals which need numerical evaluation. 
The integrands of these latter integrals are expressed in terms of recursive routines. 

The numerical integrations aie in terms of the polar angle of the wave propagation vector 
with respect to the earth's magnetic field. Singularities of the integrand are encountered when 
the applied frequency is below the plasma or gyrofrequencies. It is then necessary to carry 
out the numerical integration in the complex plane. This has not been investigated in this 
paper. Consequently the present numerical integration procedure is restricted to cases where 
this situation does not occur. The closed-form terms, however, are valid in all parameter 
ranges. 

Numerical calculations of impedance for selected values of the parameters are exhibited in 
a set of curves. 

It is pointed out that the problem may be recast variationally so that the current distribu- 
tion is computed rather than assumed. This approach, which requires no substantial changes 
in the present analytical and numerical procedures, will be given in a later paper. 
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Evaluate 

The St ai e roots of 
where 



1 1 . Appendix A 

^"S (S?-SJ)(S?-S2)0S?-S?)" (A1) 

S s i -<I 1 S'i+d 2 St-d a S 2 t +ch=0, (A2) 

d^St+Sl+Sl+Sl 

J 02 02 I 02 02 J_ 02 02! 02 02 I 02 02 I 02 02 

Cf 2 0l02T0i03T*Ji04T"0203T02^4 I & 3& 4 

,1 Q2 Q2 Q2 I Q2 Ci2 Q2 J_ Q2 Q2 Q2 

(1 3 A3 1 2 03-]- O 1 020 4 -j- O2O3O4 

,1 Q2 Q2 02 Q2 

(I4 1 02^3^4« 
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Consider the determinant 
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Q(> Ofi Q(5 or, 


4 


Q4 04 C*4 04 
0x020304 


2 


Q2 Q2 Q2 Q2 
O 1 02O304 





QO QO 00 QO 
O1AJ2O3O4 



= (S?-S|)(tf-$)(S?-SJ)-(S1-,SI)(S3-S;)(SS-,S2). 



(A3) 



This is an alternating function [Aitken, 1939] of the St, i=l, 

Let 

2n 
4 
2 




Sl n STSt n ST 


sn 


St 


Si 


o 4 


si 


si 


Si 


o 4 


si 


si 


s° 3 


04 



(A4) 



The four terms of A 2n in (Al) may be gathered over a common denominator which may be 
written in the form (A3). When this is done, it is seen that the resulting numerator is simply 
an expansion of the determinant in (A4). Thus 



2n 

4 

2 " 




(A5) 



For 27i=0, 2, or 4, two rows of the determinant in the numerator are identical and the 
determinant vanishes, so that 

A<>=A 2 =A 4 =0. (A6) 



For 2n = 6, numerator and denominator are identical: 

-48=1. 

Rearranging (A2) and multiplying through by S 2 i n ~ s , w>4, we have 

By (Al) or (A5) this implies 

-<4.2n = A.2 n = diA2 n -2 U2^2« - 4 H~ ^3^-2^-6 ^4^.2n-8j W>4. 



(A7) 



(A8) 



Equations (A6), (A7), (A8) provide the basis for recursive evaluation of the A 2n encountered 
in the text. 



Evaluation of 



12. Appendix B 
(a)/,=^J |a+6sinj8|»dl8 

i r 2ir 

( h)(Jn= 2Tj ™ s2 £|a+&sin/W. 
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Define 



fu fv fy 

I n (x,y\\a,b) = I (a+b sin P) n dl3=(i (a+6 sin I3) n - l d(3+b sin &(a+b sin f$) n - l d$ 

=al n „ l -b cos $(a+b sin &) n ~ l \ v x + 6 2 (n-l) (cos 2 P(a+b sin 0) n - 2 <//3 

2 . 2 [( a+6 sin/3)— a] 2 
cos /S=l — snr tf=l — tt 

Hence 

In=al n -i-b cos /3(a+6 sin (3) n -^ + (n-l)(b 2 I^ 2 -I„+2aI tl ^-a 2 I fl - 2 ) 

=- {(2n-l)aI n _ l + ( : rb-l)(b 2 -a 2 )I n ^ 2 -b cos (a+6 sin ^) n " 1 ||}. (Bl) 



(a) 



i f2*r i r*27r _ 

fn=2^\ |a+6sin/3| n d/3— J |a + 6 sin 0|»dj8 



where a=\a\, 6 = |6|. 

Case (i): If a>b, then a+6 sin |8 does not change sign in (0, 2ic), so that 

6 cos /3(5+6 sin p)*- 1 ! 2 ,* =0. 
Hence from (Bl), 

/„=^[(2n-l)a/ B _ 1 -(n-l)(a 2 -6 2 )/„- 2 !- 

At/ 

Case (ii) : If 6>5 ? put (B2) 

c=sin =• 
b 

Then 
2tt/„= I 2 " |5 + 6 sin 0|»(#= P (5+6 sin 0) n dp+ \ * |5+6 sin /3|»<# 

J Jo J IT 

= ( * (5+6 sin P) n dp+ \ * |5— 6 sin $\ n df$ 
Jo Jo 

- | (5+6 sin /3)*rf0+ | (5-6 sin 0) n dp+ \ * (-5+6 sin 0)*dj3+ | (5-6 sin p) n dp 

Jo Jo Jc Jtt-c 

= 2 | (5+6 sin (S) n dfi+2 f (5—6 sin p) n dp+2 | ° (—5+6 sin /3) V0 

Jo Jo J c 

=2/ B (o, | ||S, &)+2/„(0, e||a, -b)+2I n (c, w/2\\-a, b) (B3) 

^2x(L7„+F K ) 
where 

2xC7„=2/ H (o, 1 1 \a, b\+2I n (0, e\\-a, b) 



2 v V n =2I 



n (c, | H-5,6)- (B4) 
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Applying (Bl) to the right-hand side of (B3), we have for the last term of (Bl) 
-b cos p(a+b sin 0)*- 1 \l+2b cos 0(a— b sin fi) n ~ l \ c -2b cos /3(— 7i+b sin jS)*" 1 |J /2 sin p) n - 1 \V 2 
= 2ba n ~ l —2ba n ~ l +0 

Hence from (B3), (B4), and (Bl) we obtain the recursion formula for jfo in case (ii): 

f n ^l{(2n-l)aU n ^+(n-l)(b 2 -a 2 )Un-2}-l; {(2n-l)W n ^-(n-l)(b 2 -a 2 )V n - 2 } 

Ti lb 

=\ {(2n-l)a(C7 K _ 1 -F B _ 1 ) + (n-l)(6 2 -a 2 )./ n _ 2 }. (B5) 

7b 



Also 



U n =l[(n-\){V-a?)U n -,+ {'2n-\)W n ^], 

lb 



F n =i[(n-l)(6 2 -a 2 )F B . 2 -(2n-l)5F )l _ 1 ] 

lb 



Z IT 

w 

V 1 =-aV +-(b 2 -a i )H- (Be) 

7T 



Consider 



(b) g n =ir \ cos2 ^\ a + b sin P\ n <lP=7r ]* cos 2 Pla+bYm /3\ n dl3. 

&R Jo - 71 " Jo 

J n (x, y\\a, b) = cos 2 p(a+b sin p) n d&= (1— sin 2 /3)(a+6 sin 0) V0 

J x Jx 

= £ ^'-[(«+ft ft Bing)-aF j > (a+6 sin ^ 



Then we have 

Case (i) a>b: 



=p [(b 2 —a 2 )I n +2aI n+1 —I n+2 \. 



0n=p [(b 2 —a 2 )f n +2af n+1 —j n+2 ] 



Case (ii) 6>a: 

f/,=p[(i 2 -a 2 )(^+F^^ 



(Paper 68D4-355) 
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